A/ASA TTY) - 7 93 /^ 



III II 1 

III 

III 1 

III 

/b 0 

0138 4768 



NASA Technical Memorandum 79312 


NASA-TM-79312 19800005795 


MODIFICATION OF AXIAL COMPRESSOR 
STREAMLINE PROGRAM FOR ANALYSIS 
OF ENGINE TEST DATA 


Jeffrey G. Williams 
Lewis Research Center 
Cleveland, Ohio 


November 1979 


f * 

, } , » 


--...I.R 



■r 

.1 

'U 




NF00530 



E-268 


MODIFICATION OF AXIAL COMPRESSOR STREAMLINE PROGRAM 
FOR ANALYSIS OF ENGINE TEST DATA 

by Jeffrey G. Williams 
SUMMARY 

This report describes modifications of an existing 
ax ' al compressor streamline analysis computer program 
to allow input of measured radial pressure and tempera- 
ture profiles obtained from engine or cascade data. 

The proposed modifications increases the input flexi- 
bility and are accomplished without changing the com- 
puter program's input format. The computer program was 
written by Richard M. Hearsey under a grant from the 
Aerospace Research Laboratory at Wr ight-Patterson Air 
Force Base. Since this report is intended to supple- 
ment the above computer program, the reader is referred 
to Hearsey' s reports for theory, complete program list- 
ings, and detailed user's instructions. 

INTRODUCTION 

Detailed flow field information is important in 
both axial compressor analysis and design. In fan and 
compressor designs, this information permits evaluation 
of proposed improvements through efficiency calcula- 
tions. in flutter research, accurate flow information 
gives blade incidence or angle of attack which speci- 
fies the steady state pressure distribution and hence 
internal blade stresses. Likewise, in stall investiga- 
tions, accurate blade incidence is of primary concern. 

The precision required determines which two or 
three-dimensional axisymmetric analysis to use. Two 
dimensional programs (radial velocities equal to zero) 
can be reasonably accurate for hub-to-tip radius ratios 
less than 0.8. However three-dimensional analyses must 
be utilized whenever radial velocities become signifi- 
es .t. There are several such analyses that are worthy 
of consideration. Radial equilibrium is quite popular 
but requires different equations for axial stations 
located interior and exterior to blade rows (ref. 1) . 
Actuator Disk Theory, on the other hand, is rather dif- 
ficult to use in practice because of the closeness of 
'disks' in an axial compressor (ref. 2). The stream- 
line Curvature technique is extremely powerful but re- 
quires at least medium computer capability (ref. 3) . 

Of the several streamline curvature analyses in 
existence, Richard Hearsey's analysis and computer pro- 
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gram (refs. 4 & 5) can satisfy a variety of user re- 
quirements. A partial listing of these program possi- 
bilities can be found in Appendix B. The computer pro- 
gram is written in three sections which include the 
aerodynamic streamline analysis section and two blade 
design sections. All parts of the program are well 
documented and written in standard FORTRAN IV. Refer- 
ence 4 presents the theory and detailed user's instruc- 
tions while reference 5 gives program listing and two 
examples . 

Despite the versatility of Hearsey's program, the 
aerodynamic section of the program does not accept 
direct input of engine or cascade test data. This re- 
port is intended to supplement reference 4 and des- 
cribes modifications to several of the aerodynamic sec- 
tion subroutines to allow direct input of radial pres- 
sure and temperature profiles. In addition, modifica- 
tions to the CALCOMP subroutine are proposed to provide 
compatability with an earlier version of the plotting 
package. Finally, an example is prepared and run to 
sb'tw that the modifications are acceptable. 
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DISCUSSION 


Streamline Computational Technique 

Streamline curvature analysis is a means to solve 
for the velocity triangles in an axial compressor. The 
typical two-dimensional velocity triangles for rotors 
and stators (fig. 1) are extended to the radial direc- 
ti n to account for radial velocities which can be sig- 
nificant. Symmetry is assumed in the circumferential 
direction to simplify the analysis. 

The flow field grid is established by the intersec- 
tion of streamlines with specified axial stations. 

These axial stations may be inserted at strategic loca- 
tions in the flow field such as leading and trailing 
blade edges or instrumentation sites. Generally, how- 
ever, several stations must be specified upstream and 
downstream of the region of interest to define inlet 
and outlet conditions. The program will produce output 
for all defined axial stations. 

The compressible flow equations are restructured 
and written in terms of radial, circumferential and 
meridional components (ref. 4) . The meridional direc- 
tion, shown in figure 2, lies in the radial-axial 
plane. These flow equations, which contain streamline 
curvature and slope terms, are solved given an initial 
guess for the streamline pattern. The program iter- 
ates, each time improving the streamline pattern until 
the specified tolerance is achieved. Program output 
consists of velocity triangle components, flow angle, 
streamline curvature, and thermodynamic information as 
a function of radius. 

Aerodynamic Subroutine Modification 

The aerodynamic section of Richard Hearsey's com- 
puter program (ref. 4) allows a variety of input data 
cor rinations for the rotor and stator axial stations.* 
'W rk' input (DATAl) can be expressed as: Total pres- 

sure, total enthalpy, absolute angular momentum, abso- 
lute whirl velocity or relative flow angle. 'Loss' 
input (DATA 2) can be of the form: Relative total 

pressure loss coefficient, isentropic efficiency, or 


*The author assumes that the reader has an understand- 
ing of the aerodynamic (or streamline) section of the 
reference 4 computer program and is familiar witn the 
terminology and organization of the input deck. 
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entropy rise. The appropriate values of DATAl and 
DATA2 are input at each axial location as a function of 
radius or passage height. 

Typically, a researcher has correlated engine or 
cascade data in terms of radial pressure and tempera- 
ture profiles behind rotors and pressures and flow 
angles (blade metal angle plus deviation angle) behind 
stators. Figure 3 shows the typical instrumentation of 
a ''ompressor cascade. To use Hearsey's program in its 
present form the most favorable input combination would 
then be: Total pressure (DATAl, NWORK = 1) and isen- 

tropic efficiency (DATA2 , NLOSS = 2) as input behind 
rotors and flow angle (DATAl, NWORK = 5, 6, or 7) and 
relative total pressure loss coefficient ( DATA 2 , NLOSS 
= 1) as input behind stators.* This input selection 
minimizes the conversion from 'raw' data to program 
compatible data. 

Modifications can be made to the computer program 
to allow direct input of station radial pressures and 
temperatures obtained from the cascade instrumenta- 
tion. This will permit: Total pressure (DATAl) and 

to\.al temperature ( DATA2 ) as input behind rotors and 
flow angle (DATAl) and total pressure ( DATA 2 ) as input 
behind stators. Clearly these are the most advanta- 
geous input combinations when raw cascade data is 
available. The program modifications should be such 
that the original input combinations remain available 
to those researchers who might wish to use them. The 
modification scheme detailed below allows for this in- 
put flexibility. 

Aerodynamic section subroutines UD0302 and UD0307 
must be modified to allow total temperature input 
(DATA 2) behind rotors and total pressure input ( DATA2 ) 
behind stators. Subroutine UD0302 reads the input com- 
binations at each station and does appropriate units 
conversions. Subroutine UD0307 converts the station 
data into streamline entropy and enthalpy to be used by 
the main driving program, UD03AR. 

The modification concept involves sensing the mag- 
nitude of the input data and directing subroutine 
UD0307 to the correct computation of entropy and en- 
thalpy utilizing the available thermodynamic subrou- 


*The term ’relative' is misleading here but is used ex- 
tensively throughout the reference 4 program instruc- 
tions. Appendix C shows that the relative loss coeffi- 
cient reduces to the absolute total pressure loss coef- 
ficient in stator passages. 
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tines, UDGl through UDG9. For example, for stations 
behind rotors, total pressure (DATAl, NWORK = 1) is an 
acceptable input. Also, the magnitude of the isentrop- 
ic efficiency (DATA2, NLOSS = 2) should never exceed 
1.0. The modified program checks the magnitude of this 
'loss' and if greater than 1.0, assumes the input is 
absolute total temperature. Subroutine UD0307 can then 
compute enthalpy and entropy for these stations using 
th rmodynamic routines UDG2 and UDG3 respectively. 
Hence, pressure and temperature data become a suitable 
input combination behind rotors. A listing of the mod- 
ified subroutine, UD0307, which contains this limit 
checking sequence appears in Appendix D. Lines 61 - 65 
and 69 - 72 contain the specific changes for rotor sta- 
tions. 

Similar modifications are incorporated in subrou- 
tine UD0307 to handle stations behind stators. At 
these locations, flow angle (DATAl, NWORK = 5, 6, or 7) 
is an acceptable input. The relative total pressure 
loss coefficient ( DATA2 , NLOSS = 1) should be a small 
pc itive number. In fact, for efficient compressor 
blading, the relative total pressure loss coefficient 
should be less than 0.3. In this instance, however the 
magnitude of the loss coefficient does not have an up- 
per limit as did the isentropic efficiency. An arbi- 
trary limit value of 2.0 was chosen assuming the magni- 
tude of the total pressure (in any chosen set of pro- 
gram acceptable units) should never be smaller than 
this value. Again the magnitude checking procedure is 
called upon to scan for input data greater than the 
limit value of 2.0. if this condition is satisfied for 
all of the station input data, the program assumes 
total pressure as input instead of loss coefficients. 
The modified UD0307 subroutine can then calculate en- 
tropy from the total pressure input and the enthalpy of 
the previous station using subroutine UDG3 (i.e. en- 
tropy is assumed constant across stator passages). 

Sul outine UD0302, the input routine, also needs the 
li it checking sequence because a units conversion may 
be required depending on the value of the scaling fac- 
tor, SCLFAC. Subroutine UD0307, lines 142 - 146 and 
158 - 160, and subroutine UD0302, lines 190 - 198 con- 
tain the specific modifications to allow the input com- 
bination of flow angle and total pressure behind sta- 
tionary blades. A source listing of UD0302 may be 
found in Appendix E. 
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CALCOMP Subroutine Modification 

Several different versions of the well-known CAL- 
COMP plotting package exist. The main difference in 
the versions is the number of arguments in the subrou- 
tine calls to the individual CALCOMP routines. Hear- 
sey's program was written utilizing the 1970 version of 
the plotting package (ref. 6). 

The 1964-65 Version 5 CALCOMP package requires ad- 
ditional arguments in the 'CALL' statements to routines 
AXIS and LINE. The AXIS routine needs one additional 
argument, DV, which indicates the number of divisions 
per ten inches of paper to be drawn on the axis. The 
LINE routine requires four additional arguments - XMIN, 
DX, YMIN , and DY. XMIN and YMIN are minimum coordinate 
values on the graph. DX and DY represent the division 
mark incremental value for the respective axis. 

Subroutine UD0312, the plotting routine in the 
aerodynamic or streamline section of Hearsey's program, 
has several CALCOMP calls to routines AXIS and LINE. 

Tl ' s subroutine is responsible for plotting static 
pressure distributions (NPLOT = 1) and final streamline 
mesh (NPLOT = 2) . A 1964-65 Version 5 compatible list- 
ing of subroutine UD0312 may be found in Appendix F. 

Example 

Steady state data from a NASA-Lewis full scale en- 
gine test is used as input to the modified version of 
Richard Hearsey's aerodynamic program (ref. 4) detailed 
in this report and an industry developed streamline 
program. The industry program makes for a good compar- 
ison because it was developed for analysis of engine 
test data and hence accepts radial pressure and temper- 
ature distributions directly. This industry program, 
however, is not as flexible as the reference 4 program 
nor is it available to the general public. 

The analysis concentrates on the three stage fan 
mo ule at the inlet of a typical modern turbofan en- 
gine. The axial stations for computer output are 
defined and numbered in Figure 4. Stations at leading 
and trailing edges of blades follow the blade lean in 
the axial direction. Radial total pressure and temper- 
ature profiles were available behind the first two 
rotors and behind the flow straightener (F.S.). Flow 
angle was measured at only one axial location - behind 
the inlet guide vane (IGV) . The deviation angles be- 
hind the remaining stationary blades are assumed to be 
zero; that is, the flow angle equals the trailing edge 
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blade metal angle at these locations. From calibration 
plots for the flow angle probe and individual probe 
measurements, total pressure can be obtained thus al- 
lowing pressure loss input for the IGV (ref. 7). 

Figure 5 shows the non-dimensional pressure profile be- 
hind the IGV. Because of flow angle probe structural 
limitations data was only obtained in the upper 40 per- 
cent of the annulus. The remainder of the curve repre- 
sents an approximate profile. Total pressure profiles 
at other locations are shown in Figures 6, 7, and 8; 
total temperature profiles are shown in Figure 9. In 
all figures zero percent span corresponds to the fan 
hub. Approximate stator losses for stator 1, stator 2, 
and the flow straightener were computed from blade 
geometry and NACA tables (ref. 8) . 

Station input data is entered into the aerodynamic 
section of the reference 4 computer program in the fol- 
lowing manner : 

Station 1: Constant pressure. 

Pi = 136.907 KPa 
(19.856 psia) , and 
temperature, Ti = 

442. 82K (797. 08R), 
exist across engine 
inlet annulus. Whirl 
angle is assumed to be 
zero. 

Stations 2-4 No data is input 

(NDATA = 0) . These 
stations are used to 
establish the flow 
field at the inlet. 

Constant entropy 
(NWORK =0) from 
previous station (NLl 
= -1) assumed for five 
radial points (NDATA = 
5) . 

Twenty-one data points 
(NDATA = 21) specify 
flow angle (NWORK = 6) 
and total pressure 
(NLOSS= 1) as a func- 
tion of radius. 

136.392 < P < 136.795 
KPa. 


Station 5, IGV inlet: 


Station 6, IGV 
outlet : 
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No data is input. 


Twenty-one data points 
specify total pressure 
(NWORK = 1) and total 
temperature (NLOSS = 

2) as a function of 
radius. 473.82 < T < 
495.96 K. 

Notice that the modified reference 4 program is called 
upon to accept total pressure input for a stator sta- 
tion (station 6) and total temperature input for a 
rotor station (station 8). At station 6, the magnitude 
of all DATA2 (NLOSS = 1) input is greater than the ar- 
bitrary limit value of 2.0, hence total pressure input 
is assumed. At station 8, the magnitude of all DATA2 
(NLOSS = 2) input is greater than the limit value of 
l.", hence total temperature input is assumed. The 
next nine stations (9 - 17) follow the pattern estab- 
lished by stations 5-8 above. The last three sta- 
tions, 19 - 21, set the outlet conditions. Input data 
at all points is interpreted by a spline fit ( NTERP = 

0) and read in as a function of percent radial span 
( NDIMEN = 3) . A constant annulus wall boundary layer 
thickness is assumed for all axial stations. 

The subsonic solution (NMACH = 0) to the flow field 
was obtained at all specified axial stations. The mod- 
ified reference 4 computer program converged to within 
a tolerance of 0.3 percent (TOLNCE = 0.003) after 53 
iterations. Output for all twenty-one streamlines at 
each station consist of velocity components, radius of 
curvature, pressures, temperatures, and flow angles. 

The industry program was prepared and run with the same 
input information. The result was convergence to with- 
in tolerance of 0.25 percent after 38 iterations. 

Computer output for stations 7 and 8 for the modi- 
fied reference 4 program and the industry program is 
compared in figures 10 through 15. Figures 10, 11, and 
12 compare total pressure, total temperature, and flow 
angle, respectively. Figures 13, 14, and 15 compare 
the velocity triangle components: Meridional, radial, 

and relative velocity. In these figures, the modified 
reference 4 program output is symbolized with triangles 
while the industry program output appears as circles. 
All data is expressed as normalized ratios with the 
normalization constant as indicated on each graph. 


Station 7, first 
rotor inlet: 

Station 8, first 
rotor outlet: 
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Relatively good agreement is obtained between the 
modified reference 4 computer program and the industry 
program. The greatest variations in the output of the 
two streamline analyses appears at the annulus boun- 
daries. This is probably due to different methods of 
handling the boundary information or differences in the 
internal program interpolation (spline fitting) rou- 
tines. This is best evidenced in Figures 10(b) and 
l. 1 'b) . In both figures the industry program output 
deviates from both the reference 4 program output and 
the specified pressure and temperature input profiles 
at the fan hub (zero percent span). This descrepancy 
leads to differences in velocity and flow angle which 
can be seen in the remaining figures. 

The significance of this output discrepancy depends 
on the initial purpose of the analysis. As an example, 
the maximum difference in station 7 blade hub inlet 
flow angle, from Figure 12(a), is approximately 3.39° 
(-3.08 < a < 0.31°). This change in inlet flow angle 

at the hub will not greatly affect the overall blade 
pressure distribution because of the good agreement in 
flow angle output at spans greater than 10 percent. 
However, since the onset of stall flutter is sensitive 
to steady state blade incidence (ref. 9), this flow 
angle difference at the hub could be a concern.* Since 
the pressure and temperature output of the modified 
reference 4 computer program matches the pressure and 
temperature input profiles at stations 7 and 8, velo- 
city components, flow angles, and remaining thermody- 
namic output is more credible than the industry program 
output. 

Despite the differences in output from the two com- 
puter programs, a three dimensional analysis is much 
more accurate than a two-dimensional (radial velocity 
equal to zero) analysis. A two-dimensional program run 
by the author gave a station 7 flow angle of 14.15°. 

The modified reference 4 program shows that the true 
st :ion 7 flow angle ranges from -3.08° at the hub to 
14 56° near 75 percent span. The relative magnitudes 
of the velocity components, shown in Figures 13, 14, 
and 15 also indicate the significance of the radial 
velocity for this example. 


*If this difference in blade inlet flow angle had oc- 
curred at the blade tip, a more serious situation would 
result because the blade is more critically loaded in 
this area. 
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Figures 16 and 17 demonstrate the modified CALCOMP 
plotting routine, UD0312, in the aerodynamic section of 
the modified reference 4 program. Figure 16 is a 
streamline plot which shows the streamline contraction 
through the fan module. The static pressure distribu- 
tion for the hub, mid, and tip streamlines for the fan 
is shown in Figure 17. 
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SUMMARY OF RESULTS 


The aerodynamic section of Richard Hearsey's axial 
compressor streamline computer program was successfully 
modified to allow direct input of measured radial pres- 
sure and temperature profiles obtained from engine or 
cascade data. Subroutines UD0302 and UD0307 contain 
the specific modifications to permit the input combina- 
tion of flow angle and total pressure for stators and 
the input combination of total pressure and total tem- 
perature for rotors. All modifications were accom- 
plished without changing the computer program's origi- 
nal input format. 

The internal CALCOMP subroutine, UD0312, which 
plots streamline mesh and static pressure distribu- 
tions, was modified to be compatible with an old ver- 
sion (1964-65, Version 5) of the plotting package. The 
user must decide which version of UD0312 satisfies his 
computer facility requirements. 

The included example demonstrated how to construct 
th° input deck for the aerodynamic section of the modi- 
fied Hearsey program. The modified program was run and 
the results were compared with an industry program 
which accepted radial pressure and temperature profiles 
directly. Good agreement was obtained from the two 
analyses which indicates that the modifications to 
Hearsey's program were acceptable. 
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APPENDIX A 


Symbols : 

P total pressure, Pa 

Pr relative total pressure, Pa 

p static pressure, Pa 

T total temperature, R 

3 

p static density. Kg/m 

(j relative total pressure loss coefficient 

a absolute flow angle, deg 

ft relative flow angle, deg 

C absolute velocity, m/sec 

W relative velocity, m/sec 

U 'bucket' or rotor velocity, m/sec 

Subscripts : 

x axial direction 

y tangential direction 

r radial direction 

m meridional direction 

A rotor inlet 

B rotor outlet 

1-21 axial station identifier 
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APPENDIX B 


The following list highlights the flexibility built 
into Richard Hearsey's Axial Compressor Program, 
references 4 and 5. Some of the specific program capa- 
bilities are: 

1. Versatility to use program for design, off- 
design, and analysis computations. 

2. Allow intrablade station locations. 

3. Two forms of momentum equation to calculate 
subsonic and supersonic solutions of the flow 
field. 

4. Allow input in any set of consistent units. 

5. Ability to specify inlet flow distributions. 

6. Ability to run program without specifying some 
nonzero loss criterion. 

7. Multiple speed runs without inputting entire 
deck (just specify percent 'design' speed). 

8. Plots of midstreamline position showing conver- 
gence or divergence of analysis versus iter- 
ation at each station location. 

9. CALCOMP plots of final streamline mesh and 
static pressure distribution. 

10. NASTRAN compatible output (pressure difference 
for calculation of blade stresses) . 

11. Code does blade design as well as aerodynamic 
streamline calculation. 

12. Arbitrary blade shapes can be designed. 
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APPENDIX C 


In stator passages, the relative total pressure 
loss coefficient reduces to the absolute total pressure 
loss coefficient as follows: 

A subscript implies stator inlet. 

B subscript implies stator outlet. 

Symbol definitions from Hearsey's program 
instructions, reference 4 and figure 1: 

Pr ' - insentropic relative total pressure 
Pr - Relative total pressure 

P - Total pressure 

p - Static pressure 

p - Density 


pr' - Pr 



For incompressible flow: 


pr ■ = pr ft = Pft + 1/2 Pwj; 

Pr = Pr B = p B + 1/2 Pw* 

- Pr A - Pr B 

a) *= 5 — 

1/2 pw/ 

P A + 1/2 p W A 2 - (p B + 1/2 P W B 2 ) 

w = 5 

1/2 pW A 2 
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P A - P A + V2 P C a = P A + 1/2 P (CX A + Cy A } 

P fl = P A+ 1/2 p [cy 2 + w 2 - „ y 2 j 

P A - Pfl + 1/2 pfe + W A - « D - 

- 

Or: 

1/2 pW 2 = P A - p A + pU ( U / 2 - Cy A ) 
Likewise 

1/2 pw 2 « P B - p B + PU ( U / 2 - Cy B ) 
Substituting (2) into (1) : 
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For stators, u 


0, therefore: 



P - P 
A B 


1/2 P C, 


( 3 ) 


This is the definition of absolute total pressure loss coefficient. 
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APPENDIX D 


Source listing of the modified subroutine UD0307, which supplies UD03AR 
with enthalpy and entropy, appears below. Modifications are on lines 61 - 65, 
69 - 72, 142 - 146, and 158 - 160. 


SUBROUTINE U0n$07 


JO? * 


10 

11 


REAL LQSS.LAmI.LAMIPI f LA MINI SOTS 3 

COHN Ofl NST-N.S , nSTRmS.NMAX.NFORCE .NBL.NCA S C .NSPLIT ,NREAD .NPU NCH ,NPA 6$ C7 S A 

lE,NSETl f NSET2,ISTAG,ItASE,IFAlL0,iPASS, 1 ,1 YFA I L , I F F A IL , NM I X »N TR A NS* 07 S 5 

2, NPL OT ,ILO$S .LNCT.ITUB .INTO. I FA il .1 T E R , LO £ 1 .L 0 G 7 , L 0 G 3 .1 OS A .1 O s 5.l OfPTS 6 
3G6,IPRInT,NHANY,NSTPLT,NEQN 5 O 7 S 7 

COMMON NSPECl30l.NriORKlTQI.N LOS 


1.NL1 (3C),NL2 (30) ,N0IMEN( 3Q), ISl (301 , 1S2 (3 0 > , IS 3 ( 30 1 ,N£V AL ( 30 ) , ND IFS07S 

3(30 1 "noIjV^I 30 I^NCJUT iYld , N BlTaCE ( j 2 J , C U RV E C H1 3 0 * * N 

COMMON DM (11.52 ) r wFRACflI.5.2) $073 


8 

9 

10 

11 

J.Z. 


12 

13 


1A 

IS 


1* 

17 


COMMON R(21 ,301 , XL«21»30) , x< 21, 30, H (21 ,30>,S (21,301 ,VH (21 ,30),VW ($o7$ 
1? 1.3Q1 .TBETA (21 , 3 01,0 IFF ( 15. A I.POhUR (15 .A ) . FDm ID ( 1 5 . A ) . FDT IP ( 15 . A ) $075.. 
2, TER AD (5, 2) S07S 

Common oatac ( loo i .data m co , , o ata2(Ico» ,data3i 1001 .ta t aai ^ nni . D a t a s q 7 s . 
15(1001 » 0 *T a6 (100 1,QATA7( 1001 ,DATAA(1 C0> »DAtA9( 100) ,FlOM ( 10 ), SPEED (5 07$ 
23 0 ) .SPOFAC ( 10 1 ,B BLOCK ( 30) .HOIST ( 30) , WALOrK (301 .WWRL (301 ,XSTm( 15n) ,S m7$ 


13 
1A 

15 

16 
1? 

JL- 

3rstN (1501, DELF(3C),DELC (100 ),0ELTA{1CC) .TITLE ( 1 8 » , 0 ROM 2 ( 30 1 ,R IM1 ( 3$ 07$ 19 
AO ),XIHl (30) ,WORk (21 ) ,i nSS (21 ) ,t AnEPs (21 1, xl(21 1 ,VV ^ 2 1 1 , DEL >1(21 1 .L«»P7t _2 n 
5MI (21 1 ,LAMIM K21 1 .LAMIPl (21 1 ,PHI(21 1 ,CR (21) ,GA M A (21 I ,SppG(21l ,CPPG$07S 21 
6(21 ), HKFEP(21>SKEEP 121). V»KEEP (21). DELHI TO). PELT (30) $ 0.7 $ ... 2 2 


18 

19 


20 

21 


22 

23 


COMMON WiSr , SHA pE, sCLF AC, EJ.G, T 0LNCE ,XSCALE ,PS C ALE , PLO U ,R LOW , 
1 tRCONST ,FM2 ,HMIN ,C1 ,PI ,CONTR , CO N M X 


XHMAX$0?S 23 

$07$ 24 

$075 25 

$07$ 26 

S o7 $ 2 7 

*07$ ...28 

$n7$ 29 

$07$ 30 


2A 

25 


L1=I«N L 1(I) 
L 2- 1 *NL 2 ( I 1 


26 

27 


IU=NbORK(I ) 
ILzNLOSS(I) 


28 

29 


XNrSPEEO(I)*SPDFAC(ICASE }*PI/(33.0*SCLF AC » 

GO T0(1 00, 250,270, 290, MAO. A AO. A AO), IW 


30 

31 


ICO GO T0(110, 190, 210,110), IL 
1 1 0 I F(L2,NE,I)G0 TO ISO 


$07 $ 
_$ 07$ 
$07$ 
_S 07$ 
3 07-; 
tO 'S 


Jl 

32 

33 

3 1 " 


32 

33 


00 1 A 0 Jrl.NSTRMS 
IF( IPASS.E0.1. ANO. ITER.EO. 0160 TO 120 


3A 

35 


IF ( ITER.EO.O 1VV ( J) = VM ( J, I ) 

XlrH (J ,1 )-(VV ( J1»*2-»VW (J, I)»*7)/(2.0*G«EJI 


36 

37 


X2=H(J,I)-(VW(J,II**2-(VU(J t I )-XN*R( J,I )>**2)/(2.C*C*EJ) 
I F ( X 1 « LT , HMI N I X 1 -HMIN 


307$ 37 

$075 _ 38 
$07$ 39 

$075 AO 

$ 07 $ Al 
$07$ A 2 


38 

39 


IF1X2.LT.HM IN 1X2 rHMIN 

X 3 = 1 .0 / (1 .0* LOSS (Jl ♦( 1 ,0 — UOG A (X 1 ,S( J ,1) l/UDGA (X 2.S 1J.HH) 


AO 

Al 


A2 

A3 


AA 

A5 


120 
1 3 1 ' 
1A 


GO TO 130 
X3rl .0 


H (J,I >=UDG2(S(J,L1) .WORK ( Jl/X 3) 
S ( J,I)-UQG3 1 WORK (Jl ,H ( J, I 1 1 


150 


$07$ A3 

$07$.... AA 

$07$ A 5 

$ P7jt A 6 

IF(IPASS.EQ.1.AN0.L2.GT.I 1G0 TO 160 3075 A7 

Xl-Hlj,Ll>-<VW(J,Ll t**2 - ( vw( J.L 1 )~XN »R( J.LH )*w2 1/ 1 2 .0>6«E J 1 « XN*» 7$ p7$ a 8 


GO TO 230 
00 180 J- 1 .NSTRMS 


A6 

A7 


A8 

A9 


50 

51 


52 

53 


1* (R(J,II**2-R(J»L11**2)/(2,0*6*E») 
IF(X1.LT.HMIN|Xi=HMIN 


X2rH(J,L2 )-( VM( j,L2I**2*VU(J,L2)**2) /(2 .0*G»EJ) 

X3=H(J,L2>-( VU1J.L2 i*»2-(VW( J.L 2 )-XN«Rl J,L2)1 ♦»?)/( ?.0»6»EJ) 


IF ( X2.LT.hMiNIX2-HmIN 
IF(X3.LT,HMIN1X3 -HM I N 


$07$ A9 
* 0?$_ J>d 
$07$ 51 

_$ 07$ .52 

$07$ 53 

$07$ 5A 


5A 

55 

56 


X A — 1 .O-LOSS ( J)/U06A (X 1 ,S ( J.L 1 1) *(UDG tf ( X 3, S( J.L 2 ) 1 -U06A ( X2 , S ( J 

Jl 

GO TO 170 


,L2))$07$ 55 

SOTS . 56 

$07$ 57 
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5 07 




IF< iPAS 
IFIITER 
X lrH < J , 
X 2-H I J , 
IF! XI. L 
IFIX2.L 


x3=i.a/ 
GO TO 3 
X3rl .a “ 
S I J,I )r 
GO TO 5 
00 39Q 


IF( IPAS 
X 1 -H (J, 1 
1 * I R <J,I 
IF (X I. L 
X 2 rM IJ, 
X 3 -H |J , I 

I F ( X 2 • L 
IF1 X3.L 
XH=1 .0-1 


S. EQ.l.AN 
,EQ .0 >VV t 
I|-«VV(J) 
Tt-lyytJ, 

T. HM IN )X 1 
T.HMINIX 


II .C*LOSS 
6 0 

U OG 3 ( X 3 »U 
70 

Jr 1 .NSTRH 


S. EQ.l.AN 
L I )- I VM < J 
I* * 2 -R (J, 

T. HHi mXl 
L 2 I-IWMIJ 
L 2 I-IVMIJ 
T.HHINIX2 
T.HMINIX3 
LOSS I Jl/U 


s s07s S.L. 

XN*IR I J,I l*VU CJ ,1 l-RIMl (J )*VU( J ,1-1 »)/IG*EJ) 5075 82 

0.390), IL s07j 

370 5 

s s 


O.ItER.EQ.OJGO TO 35 G J0>-3 G-, 

J )- VH IJ, I ) . 0 7 5 E 7 

**2*VWIJ, I )**2) / { 2.0*G*EJ ) 5 07 S SB 

I)«»2-t'/UiJ t I )-XH»R 1 J .1 )) »*2) / 1 2 . n » G*E J ) 5 0 jrs 8 9_. 

=TMIN 507s 90 

-HH I N 5075 


(J)*(1.0-UDG4 (X l.SIJ .1) ) /UDG 14 1 X 2 . 5 IJ,I) I) ) 5 075 92 

10 7.5 9 3_ 

5075 SR 

DG4 (H ( J, I ) ,S t J, L 1 ) i ,HtJ,D) 507 5 95 

5075 96 

S 5075 97 

0.L2.GT.HG0 TO 38Q 5075 9g 

.L1»«»2-IVMIJ.L1)-XN«R( J.LH>»»2>/I2.0*6»EJ)«XN»«2S075 9 9 

L 1 )** 2 )/| 2 . 0 *G*EJ) 5075 loo 

-HH I N 5075 101 

,L2»**2«VUI j,l2 1**21/12 .q*G*EJ) 6075 102 

.t 7 )*» 2 -ivy| J,L 2 )-XN«R 1 J.L 2 »I*» 2 )/I 2 . 0 *G*EJ) 5075 103 

-HHIN 5075 104 

-hhI n s.Q 7 .a.J.Q5 . 

D 64 ( X 1 ,S ( J,L 1 )l *lU0 6 4IX3.SIJ f L2))-U0G l )lX2,St J,L2> >5073 106 


OGRIHIJ.H.SIJ.LDl.HtJ.I ) > 
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S ( J,IlruO 63 (U 0 G 4 (He J,L 1 | *L 0 SS (J >* CH { J ,1 } -H C J, L 1 ) ) , S ( J,L 1 1 > ,H C J, I 11 5 075 

GO TO 570 _ _ S. 07 S 

DO 43 o J= 1 ,NSTRMS 5075 

S ( J,l> =S f J.H|«L OSS ( J 1 5 Q 7 5 

GO TO 570 5071 

0 ii 50 J =1 .NSTRMS 






GO TO <480 

G X 3-1 . 0 

0 S ( j tXJ-U0G3 (X3*UCG4 ( X 2 f S ( J »L 1 1 J ,X2> 

GO TO S 70 

*490 K PRES rO 

DO <4 52 J=1 .NSTRH 


492 IFUBS(L0SS(J)).GT.2.0»KPRES = kPP£S*1 

IF (KPRES.EU.NSTRMSIGO TO ; d; 

DO 5 00 J=1.NSTRHS 

X4 = X1( J >♦ ( XN»R ( j,H )*»2/(2. Q «G»EJ> 

IF 1 X4 .LT .HHIN IX 4 XHMIN 

ylrUQG4(X4.S( J.L 1 >> 

If(IP*SS.EQ.l.AND.L2.GT*I IGO TO 500 

X2=XI(J)*(XN»R( J.L2I >»*2/(2.0« G*EJl 

X3zH ( J,L2 )- ( VM( J ,L2 l**2» VU { J,L2 »**2» /(2 .D*G*£J I 

IF ( X 2.LT.HMIMX2-HMIN 

IF(X3.LT.HMIS)X3-HHIN ~ 


5n75 130 
_5 0 7 5 131 
5075 132 
_5 07 5 133 
MOD . -JGW 
HOD . -JGW 


MOD . -JGW 
.MOD . -JGW 
MOD . -JGW 
_S 07 $ 135 
5 075 136 

5 0 ~ i •- 

C fj 7 *, 5 •* £ 

.5075 13R 
5075 140 
.5.075 1 4 1 
5075 142 
5075 143 


a 500 

S ( J , I ) rUDG 3 ( X 1 t X 4 ) 

5075 144 

T 

GO TO 570 

507$ 145 

5 05 

DO G 06 J=1,NSTRMS 

MOD . -JGW 

) >06 

S ( J „l I-UDG3 (LOSS ( J) ,XI (J J ) 

MOD . -JGW 

3 

GO To 570 

MOD .-JGW 

l &1Q 

IF(IP*SS.E0.1.AND.tTER.EO.0»GO To 530 

5075 146 



DO 520 J- 1 t NSTRHS 5075 147 

IF(ITER.E0.01VV(J| = V Mj J.Il 5 075 14 8 

Xl-HCJ.I-l )*XN*( VV( JI*(TBE tA( J,I»«xN*R( J.IJ/VV (J) )*R(J,I)-RIM1(J)*!075 149 

1WWI J.I-l l 1 / ( G<E J 1 * 07 5 150 

IF(X1.CT.HMINIX1:HMIN 5075 151 
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Source listing of modified subroutine UD0302, the input routine, appears 
below. Modifications are on lines 190 - 198. Statement labels which indicate 
a date (line 177: 11/30/77) are program updates. 


s ubroutine: U 00302 




REAL LOSS,LAMI,LA*IPi,LANIMl J P 25 

COMMON NSTNS T NS TRMS .NMAX , NEORCE . NE L.NCAS E t N S ?t I T .NR EAO . NPUNr H . NP A &S 0 2 5 
IE ,NSET1 ,nSET2,I S TAG , ICASE ,IF A ILO.IPA SS, I, IVFAIL .IFF AIL,NMIX,NTRANSS025 
2,NPL0T ,ILOSS ,LNCT,ITU9,IMI0. I FA IL»ITER«LOG1,L062.L003,L064«LOG5»L05025 


10 

11 


12 
1 3_ 
L*» 
15 


16 

17 


3G6, IPRlNT.NMANy.NSTPLT ,4EQN 5025 7 

COMMON NSPeC< 3Q),NMOR» (( 3Q>,N L('SS(3G1 ,NOAy A_! 3j)J , NT ER P ( 30 I , N H A CH I 30 >5 02 «. 8 

1.NL1 13G».NL2<3G) » NO I HEM JO) , I S ; ?C ;; :• 2 { 3 01 , jS 3 < 30 ) , N E V AL (3o> .NDIFS02S 9 

2F <91 , NOEL < 30 l.NL ITER < 50 ,NM t 2 J ^ i « D( ? ).NCURVEC 3Q.NWH ICH 1 3 0 > . NOUT I * 0 2*_ Jl 0 
3<30),NOUT2(3C),noUT3(3g),NBLADE(30) %02 4 11 

COMMON 0M(11,5,2> , U FR.M11.S.2I 50 25 1 7 


COMMON R(21,3G),XL(21,3g),XT21* 30 > , H I 21 , 3 0 ) ,S ( 2 1 , 3 0 ) , VM ( 2 1 , 30 > , V W ( 5 025 
121, 3-1 , T BEtA( 21,301,OIFF115.8 ) r FQ H UB (15.41.F0MID11 5 .4) .FOTIP ( 15, 4 1 % 02 % 
2. TER AO (5,21 5025 

COMMON OATACdnnl.OATAKlQul.OATA 2(1 C C 1 ,nATA 3(1001 , D AJ A9J 1 001 ,0ATA5 P 25 
15(l0C,l,OATA61100),0A7A7(10C),0ATA8!lC0l,OATA9(iO0l .FLOW ( 1 0 > , SPEED ( 5 02 5 
230 ) ,SPDFAC( 101 .BBlOCX ( 30 > ,BOIST 1301 , wBLOCK( 30) ,wyBL (301 ,XSTN( 1501 ,^02$ 


13 

14 

15 

16 

17 

18 


18 

19 

20 


3RSTN I 150) ,DELF< 3CI,DELC< 100) , DELTA! 1 00) ,TITLE( 1 8 1 ,DR0M2 ( 30 1 ,R IM1 ( 35 02 5 
4Q 1 , X IM 1 (30) , h OR 8 (21 ) .LOSS ( ?] 1 ,T ANEPS ( 21 ) , XI ( 2 1 ) , VV ( 21 ) , 0FLM12 1 )_, l» * 0 2 5 
5 M I ( 2 1 ) ,LAMIMl(21 l.LAMIPl (21) ,PHT (21 ) ,CR (21),GAMA(21 ),SPPG(21) .CPP65025 


19 
.2 0 
21 


22 

23 

COMMON VISK, SHAPE, SCLFAC,EJ,G ,T OLNCE ,XS C* LE , PS C AL E , PLOW 
1, rCONST ,FM2,HMIN ,C1 ,PI ,CONTR .CONMX 

,RlOW,XmMAX5o25 

5025 

23 

29 

24 

DIMENSION I j (21 ,30) ,JJ(21,3o) 

5 02 t 

25 

25 

ECUIVALENCE ( H , 1 1 ) , ( S , JJ ) 

5025 

26 

26 

COMMQN/PAGE/LIMI T,LO 

S 0 25 

27 

27 

NEVAL( 1 >=0 

5025 

28 

28 

RE*0 (LCG1 ,10C )T iTpE 

502$ 

2° 

29 

100 FORMAT (13A4) 

5025 

30 


5025 30 

5 025 31 

DATA,/, IPX, 10(1H*),// ,10X,5HTITLE,34X,2H- .18A5025 32 

502* 33 

5025 34 


30 

31 
32' 

3 3_ 

39 

35 


110 


14 


LNCTrLNCT ,4 

r * i i .inrtfi. 


CALL UOGl(LNCT) 5025 

REA0(L0G1,12G>NSTNS.NSTRMS,NMAX .NFORCE, NBL.NCASE ,N S PL I T ,N S E T 1 ,NS E T5 0? 5 


35 


36 

37 


38 

39 


40 

41 


120 


12.NREAD ,NPUNCH,NPLOT ,NP A G E ,N T R AN S ,NM IX, NmA*Y , NSTPLT ,NEQN 
FORMAT) 1813) 


IF(NSTRmS,E0»C)NSTRMS-11 
IF (NMAX .EQ .0 INM AX=40 


SC2‘. 

■025 38 

5025 39 

5025 40 


IF (NFORCE. EO .n) NFORCE =10 
IF(NCASE.E0.C)NCASE=1 


5025 

5025 


4 1 
42 


42 

43 
"44 

45 


46 

47 


IF(NP*GE.E0.P)NPAGE^80 
L Q = LGG 2 


5025 

5025 


43 

44 


L IM I T-NPA GE 

CALL U00303(LNCT,19) 


5025 45 

5025 46 

WRI TE1L0G2, 1 30INSTNS.NSTRMS , NMAX, NFORCE ,NBL,NC ASE ,N SPLI T, NSET 1 ,NSE5 025 4 7 

1T2.NREA0.NPUNCH.NPLQT , NP A GE , N TR «NS, N MIX ,N MANY , N STPL T ,NEQN 5025 48 


48 

49 


50 

Jl 

52' 

53 


130 


FORMAT (2X,/, 10X, 1 &HNUH8E R OF ST A T ION S ,2 IX , 1 H= , I 3 , / , 1 OX , 21 HNUMBER 05025 49 

IF STREAMLINES , 1 8 X ,1H = , 1 3 , / , 1 0 X, 20HM A X NUMBER OF PASSES, I9X . 1H =i_I 3 j* 0 2 5 50 

2/.10X.30HMAX number OF ARBITRARY PA S SES ,9 X, 1H = , 13 , / , 10X ,2 9HB0UND AR* 02 5 5l 
3T LAYER CALC INDICATO R , 1 0 X, 1 H =, 1 3 ,/ , IPX ,2 4HNU W B ER OF RUNNING POINT* 02 5_ 52 
4S,15 x, 1H = ,I3,/,10X,33HSTREAMLINE DISTRIBUTION INDICATOR , 6X,1H = , I 3, i 025 53 

5/ , 1 OX, 34H<|UmBER OF L 0 s S/ 0 -FA C TO R CURVE SE TS ,5 X . lHr , 1 3 ,/ ,1 0 X , 34 HNUMj 0? 5 54 


54 

55 


56 


6BER OF LOSS/ T.E . LOSS CURVE S E TS ,5X, 1 H=, 13 ,/ , I □ X , 26H STRE AML INE INPUSC25 55 
7T INDICATOR ,13X,lHr, 13, /.IPX, 27HSTREAMLINE OUTPUT INDICATOR , 12X.1H5025 56 
8-,I3,/,10X,24HPRECISI0N PLOT IN 0 ICAT OR , 15 X , 1H= , 13 , / , 10X ,2 9 HM A X NUM5025 57 



57 


9BER OF LINES/PAGE, 15 X .1H=, 13, /, 10X,29HuAKC TRANSPORT CALC 1NDICATOJ02* 

58 

58 




jJS 025 

. 59 

59 


B,10X,33HN0 OF STATIONS FROM ANALYTIC SECfJ, 6 X, 


60 




61 

61 


D1NDICAT0R,7X,1H=,I3> 


5025 

62 

62 


ITUBzNSTRMS -1 


»02 5 

63 

63 


I HID -NS TR MS/ 2*1 


502* 

69 

69 


iFlNMANV.LQ.OGn To 13b 


Sf)25 

65 

66 

65 


READlL06l,132)CNNHICHm ,IZ1,N« '5 


*025 

66 

132 

FORMAT 129131 


1025 

67 

67 


CALL UD03031LNCT ,2) 


502* 

68 

68 


«riteilog2,139m nvhici: : r » ,1 = 1 .nhany i 


5 fl?J 

69 

69 

70 

139 

F0RMATljX,/,10X,51HGE0KETRY c° m ES from analytic SEC jION for STATI0S02S 70 
INS ,2313) 5029 71 

71 

7? 

136 

CALL UDp3 0 ^ ( LNC t , 71 

REA0IL0Gl f 190lr.,FJ.SC|FAC,TnlfjCE,VI^K,9HAPE 


5C2* 

5025 

72 

73 

73 

190 

FORMAT (6F12 , c ) 


*025 

79 

79 


IF1G.E0.0. DIG-32 .179 



75 

75 

76 


lF<EJ*EO.O.O)Ej-778.16 

IF(SCLFAC.EQ,0.0)SCLFAC-12.0 


*C 2 * 

502* 

76 

77 

77 


IFtTOLNCE.EQ.Q.O ITOLNCEzO.OOI 


5025 

78 

78 


TFIVISK.Eo.n.nl VISK- 0.000 18 


5_02_5- 

79 

79 

80 


IF«5HAPE.E0.C.0)SHAPE-0.7 

WRITECL0G2, 1501G ,EJ,SCLFAC,TOLNCE ,VI SK, SHaPE 


5C25 

5025 

80 

81 

81 

82 

150 

FORMAT (2 X ,/ » 1CX t 22HGRAVI TATIONAL CONSTANT ,17X 
1ULES EQUIVALENT , 22X, 1H-,E8.3 ,/, 1 oX,29HL1NEAR 

, 1Hz,F8.i,,/,10X,17HJOsO2S 
DIMENSION SCALE FACT0SC2* 

82 

83 

83 

69 


2R,10X,1H=,F8 .9,/ , 10 X, 1 BHBASIC TOLERaNcE ,2 9X , 1 H = ,F8 . 5 ,/ , 1 OX , 1 9 HK INES 025 
3MATIC VISCOSITY, 20X,lHz.F8. 5. /, 1GX, 1 7HB .L . SHAPE F A CTOR , 22 X, 1 H= , F 85 02 5 

69 

85 

85 

86 


9.5) 

CALL UDQ3031LNCT ,71 


*025 

5025 

86 

87 

87 


READ(L0G1 ,19 0 IX SCALE , PSC ALE , RLO V , PL 0 V ,X MM A X ,R CONST 

* 02 * 

88 

88 


IFCXMMAX.EU.C.O) XMMAX-O . 6 


5025 

89 

89 


IF CR CONST .EQ.0.0 1PC0NSTZ6.0 


SC25 

90 

90 


URXTE(LOG2,16q)XSCALE,PSCALE,RLOU,PLOV, XMMAX, 

R CONST 

5 02* 

91 

91 

160 

F0RMAH2X,/,10X,29HPL0TTING SCALE FOR 0 IMENSI ONS , 1 0 X, lHr, F 7. 3 ,/ , 

10502* 

9? 

92 


IX »28 H PL0TTImG SCALE FOR P RES S UR FS , 1 1 X . 1 H r ,F7. 

3 ./.10X.22HHINIHUM 

R A S 0 ' 

s • 

93 


2DIUS ON PL0T,17X,1HZ,F7.3,/,10X,29HMINIMUM PRESSURE ON PL0T.15X, 

1HI026 

99 

99 


3-,F7.3,/,l DX .9 OH MAX I MUM M-SQUAREO IN RELAXATION FACTOR Z.F8.9,/, 

105025 

95 

95 


9X ,29HC0nSTANT IN RELAXATION F AC TOR, 1 OX, 1Hz ,FB 

.9 1 

5025 

96 

96 


CALL UD0303 (LNCT ,3) 


5025 

97 

97 


READ ILCG1 ,190)C0NTR,C0NMX 


5025 

98 

98 


WRITE (L 062,1 69 > CONT R , CON MX 


5021 

99 

99 

169 

format»2x,/,iox,22hwake transfer constant, 17x 

, 1HZ.F8.5 ,/, 1 CX, 25 HT Us 021 

ICO 

ICO 


1 RBULENT MIXING CONS T AN T, 1 9X , 1 Hz , F 8 . 5 1 


* 02 * 

lei 

101 


CALL UD0303 (LNCT ,5*NCASE 1 


5021 

102 

102 


READILOGI ,170 (FLOW (K ) .SPDFAC (K 1 ,Kr 1 ,NCASE) 


5025 

103 

103 

1 0 

FORMAT I2F12 .01 


5025 

109 

109 


WRITEtLOG 2 , 180 ) (K,FLOU(K 1 ,SPDFaC<K) ,KZ1,ncaSE) 

5025 

105 

lcs 

106 

180 

F0RMATI2X,/ ,10X,21HPOINTS TO BE COMP UTE 0, // , 1 0 X ,2HN0 , 6 X , 8HFL0W R A T E* 02 5 
1,9X»12hSPEED FACTOR, //,C 1 Ox . I 2. F 1 3. 3 ,F1 9. 3 ) ) S02* 

106 

107 

107 


READ 1L0G1.19C1L1 ,CXSTN(K ) .RSTNIK) ,K = l,Lll 


5025 

108 

108 

190 

FORMAT (13,/ , I2F12.0) > 


1025 

1Q9 

109 

110 


I S T A 6 — 0 

IFCRSTNU).fQ.0.0)ISTAG=1 


* 02 * 
% 02 S 

110 

111 

111 


N SPEC 1 1 1 zLl 


1025 

112 

112 


CALL UD03031LNCT »7*L 1 1 


5025 

113 

113 


UR1Te«L0g2,2DQ) l1»<*STn«K),RsTN IK),K=1,L1 1 


I02S 

11 “ 


22 


114 

200 

FORM*! 12X,/,10X,36HAknULUS / COMPUTING 

STATION GEOMETRY 

.//,10X,24Ht02j 

115 

115 


15 T A T ION 1 SPECIFIED BV.I3.7H POIN T S ,/ /, 17X. 4 hXSTn ,8X . 


116 

116 


2F22.4.F12.4) ) 


5 02 1 

117 

117 


IS1 111=1 


SO?* 

lift 

118 


L AS T =L 1 


S 02 S 

119 

119 


DO 220 I=?,NSTNS 


S 021 

120 

120 


READ (L0G1 ,21C»L1 


S02S 

121 

121 

210 

FORMAT <131 


so?s 

.122 

122 


NEXT:LAST*1 


s 02 s 

123 

123 


LAST=LAST*L1 


sn2s 

124 

124 


IF<LAST.GT.150»G0 TO 550 


S 02* 

125 

125 


REaoILOgI *170) IXSTfj : .oSTNCK 1, K=N£xT,L ASTI 

*02* 

l26 

126 


IF(RSTN<NEXT >.EO.O.C)ISTAG=I 


*02* 

127 

127 


CALL UD0303|LNCT .5«L11 


*021 

128 

128 


IS1 <I)=NEXt 


*021 

129 

1?9 


NSPEC1I 1 =L 1 


*021 

130 

130 

220 

WRITE <LOG2,2 3o) I ,Ll , <XSTN <K) ,RSTN <K ) ,K = 

NrXT.LASTI 

*02* 

131 

131 

23q 

F0RMATC2X,/, 10X.7HSTATI0N.I3. 14H SPECIFIED BY.Ix.tH PO IN TS , / / . 1 7X* 02 * 

132 

132 


l,4HXSTN,8X,4HRST Nt //, <F22.4,F12 .41) 


*0 2 * 

133 

133 


SPEED! 1 )=D.O 


*02* 

134 

134 


READ! L0G1 , 240 )L1 ,NTERP<1 ) ,NDI ME N < 1 ) , NM . CH ! 1 ) , < 0 AT AC <K » • 

DATAltK),DA*o2* 

135 

135 


1TA2(K),DAtA3CK),K=1,L1 1 


*021 

136 

136 

2 ^ 

FORMaT <413,/ , C4F12.0 ) ) 


*02* 

137 

137 


CALL UD03C3<LNCT,7«L1 ) 


*02* 

138 

138 


I S2 < 1 ) = 1 


502* 

139 

139 


N CA T A < 1 ) = Ll 


*02* 

IMP 

140 


LAS T=L 1 


*02* 

141 

141 


WRITE<L0G2,2SD)L1,NTERP< 1),NDIMEN<1 1 .NMaCHU), <DATAC<I<I 

,0 AT A 1 < K ) ,0* 02* 

142 

142 


1ATA2!K),DATA3<K),K = 1,L1) 


*02* 

143 

143 

250 

F0RMaH 2X,/ ,10X,24HST»Tl0N calculation 

DATA,// ,7X,18HSTATI0N l NS02* 

144 

144 


10ATAr,I3,7H NTERPr f l 2 , 8 H ND IMEN = , 12 , 7H 

NMACH=, x2,// »11X ,5HDATAC.6X*02* 

145 

145 


2, 14H TOT AL PRESSURE, 4X.17HT0TAL TEMPE R AT UR E . MX , 1 lHUHIRL 

ANGLE,//, <5*021 

146 

1 4 6 


3X,F12.4,f15.4,F19.3,F18.31) 


* 02 j 

147 

147 


00 252 «=l,Ll 


*02* 

148 

148 

252 

DATAKX )zDATAl|K )*SCLFAC**2 


*021 

149 

149 


LAS TD=0 


1021 

! r ' 

150 


NOUT 1 1 1 ) =0 


* 02 * 

15; 

151 


N0UT2 < 1 1=0 


-02* 

152 

152 


DO 320 1=2 ,NSTNS 


*02* 

153 

153 


LOgN=LOG1 


*02* 

154 

154 


IF <NMANY.EU.O)GO TO 258 


5021 

155 

155 


DO 254 Llrl.NMANY 


*02* 

156 

156 


IF<NwHICH<Ll ) .EQ .1 )G0 TO 256 


102* 

157 

157 

2T' 

CONTINUE 


*02* 

15.8 . 

158 


GO TO z58 


102* 

159 

159 

2 - 

L0GN=LCG5 


1021 

160 

160 

2 J 8 

READ <L0GN,260)NOATA 11) ,NTERP< I) ,NDIMEN< I ),NMACH (I ) ,N»ORMI > ,NL0SS<*02* 

161 

161 


lD.NLKl) ,NL2<I 1 ,NEVaL<I ) ,NCURVE< I) »NU TER < 1 1 , NOEL ! 1 1 * N OUT 1 CII , NOU* 02* 

162 

162 


2T2 < I 1.N0UT3 1 I ) ,N BLADE < I) 


*02* 

163 

163 

260 

FORMAT <16131 


1021 

164 

164 


Ll = 3 


5025 

165 

165 


IF(N0ATA<I|.NE.0)Ll=Ll*5*NDATA<I) 


*021 

166 

166 


IF 1MDEL <I).NE.0)LlrLl«3*NDEL < I ) 


*021 

167 

167 


Call uq0303 ilnct ,l1 ) 


*02* 

168 

16 6 


WRIVE1L0G2 ,27 0)1 ,NDATA <1 ) .NTERP <1 I.NqIM EN < 1 1 , N H ACH < I ) , NW 0 RK < 1 1 , NLO* 02S 

169 

169 


ISS < I ) ,NL1 <1 ) .NL2 <1 ) . NEVAL <1 ) . NCURVE < I ), 

NLITER < I I.NOEL! I ) , NOUT U 1 1 , S 02 * 

170 

170 


2N0UT2<1),N0UT3<I >,NBLADE<I> 


5 02S 

171 


23 





270 F0RMATI2X,/, 7X,7HSTATI0N,I3,8H ND A TA= ,1 3,7H NTERP=,I2,8H N0IKEN=*p2S j72 
1.I2.7H NMACHr.l2.7H NU0RKr.T7.7H NL0SSr.t?.5P ».,lr.T3.SH Nl 2= . T 1 . 7* n? 5 I7T 


2H NE VALr, 12 f 8 H NCURVE=,I2 , 8 H NL I TER r ,1 3 , fa H NDE L = , I 3 ./ , 1 9X , 6 HN 0 UT 1=S02* 171 


SPEEOI 1 1-0.0 5025 176 

IFiNOATAlII.EO nlGO TO 720 . _ * 


IFINUORKII) .E0.7INL0SS jl jrl 11/30/77 

NCXT-L AST*! 102s_1M 

LASTrLAST«NOATA I II *02* 179 

IS2ll)=NEXT *02* IB q 

IFILAST.GT.IOOIGO TO 550 *02* 181 

REAO<L0GN,2aC)SPEEimi.<DATAC IK S .OAT AUK) .DATA 


IK 1,0 AT A5 IK I, OAT A 6 IX 1 , 0 AT A 7 1 K » , D A T AS I K I , PA T A9 1 K » ,K=NEXT , LA ST » *02* 183 

FOR MAT I FI 2.0,/ , I6F12.0./.9F12.0II * 02* 1 8* 

UR 3 TE 1 L 0 G 2 .2 90 IS PEED II),(0ArACIKI,DS.TAllK|,0ATA2tKJ .DAT A3 IK I , OAT A95 02* 185 
llK| .OAIAsIK) , 0 A T A 6 I K ) . D» T A 7 I K I » n A TA 8 lnl . 

FORMAT (2X,/,10X,7HSPEED r,F9. 2, // , 1 3 X . 5 HD ATAC * 7 X, SHOaT a 1 1 7X ,5 H 0 AT AS02* 187 
12 , 7 X, 5 HO AT A 3 ,7X , SHD A 7 A 9 , 7 x , 5 HDA TA 5 , 7 X ,5 HD AT A 


28 .7X.5HDATA9.// . I 10 X , F 9. 9 , f 1 2 .3 , F 13 . 6 ,F 11 . 9 ,F1 2 .5 , F 12 .5 ,9 F 12 .9 1 > *D2* 189 

IFiNLOSSlll.NE.llGO TO 29 


IFIKPReSS.NE.kCHECkIGO To 293 
DO 292 K rNEX T .LAST 


292 DATA2lKlr0ATA2(X )*SCLFAC**2 

293 IFIN-OK* 111 . n E.1 IGO TO 296 


DO 29i, K=NEXT,LA$T 
>9 DATA1 (K irOATAlIK )*$CLFAC**2 


16 IFInEVALUI.GT.O.AND.NSTRMS.GT. NDATA 

IFINDEHII.rC.0lGO TO 320 

NEXT rL A ST D * 1 
LASTOrLASTO.NOEL III 


I S3 1 1 ) rNEXT 

IFILASTD.6T. 100160 TO 550 


RE AO IL 0 G 1 . 300 U 0 ELC IK I .DELTA I K| ,KrNE 
}C FORMAT I2F12.C1 


URXTE|LOG2,31ol I DELC |K | , DELT A |K I ,KrN 
tQ FORHAT(?Xf/.13X.9H0ELC,8X .5HDELTA,// 
?G CONTINUE 

CALL U 0 c 3 C 3 ll.NCT ,S»NSTNS 1 


READ (L0G1 ,330 UNBLOCK |1),B8L0CK|I),B 

FORMAT I 3F12 .01 

URITEIL062,39Q) 1 1 ,W BLOCK 1 1 1 , B BL OCX I 1 
0 FORMAT ( 2X ,/ , 10X , 3CHBL0CK A GE FACTOR S 


IN NALL BLOCKAGE uAKE BLOCKAGE 
23X,m,F16.5,F16. 5.F19.31 1 


XFINSET1.EO.OI60 TO 380 
DO 370 K=1,NSET1 


REA0lL0Gl,350lLl.lDIFFlJ,K),FDHUBIJ, 

U 

FORMAT 113,/, (9F12.0I I 
ALL UD0303ILNCT »6«L1I 


URITEIL0G2.360IK ,L1 , ID IFF I J,K », FDHUB 

11 . Ll| 

>0 format i 2 x,/,iox.55hloss parameter / 


M 0 D.-J 6 U 

MOD.-JGU 


mod.-jGU 

MOD.-JGU 


*025 191 
*02* 192 


III )LAST=LAST«NSTRMS-NOATAIXI *02* 193 

502* 1 99 

*02* 195 
*n2* 196 


*02* 197 

*' * 


XT.LASTDI 502' 1 y9 

iC2S POD 


EXT.LASTd* * 02 * 201 

,(1CX,F9.9,F12.91> 50 25 202 

*025 2o3 
5025 209 


Fl2. 91) 


DISTII),I=l,NSTNS) 5025 205 

1 0 2 * 20 6 _ 

"l.BOISTIII ,I = 1,NSTNS1 *02* 207 


I.BDISTIl) .I-l.NSTNSI 502* 207 

PEC IF1CATI0NS,//,10X,66HSTATI050Z* 208 


WAKE DISTRIBUTION FACTOR,//, 1 1502* 209 

502* 210 _ 


5025 211 
5025 212 


KI,FDMIDjJ,K) f FOTIP (J,K», J=1 ,L5C2* 213 

5 JO 25 219 

5025 215 
5025 


I J.Kl ,FDHI0IJ,K l.FDTIPC J,K», J=5025 217 

50 25 218 

DIFFUSION FACTOR CURVES FOR BLS02« 2X9 





228 


1 A DC TYPE » 12 , 1 5 » 1 6H D -FACTORS 61 VEN,/ / ,1 5X ,9HDI FFUS 1 ON ,5X , 3CHL 

0 S *02* 

220 

2?9 


25 Parameter s,/,16X,7hfactop* ,g~.3hhub,9x ,3Hmid.8x 

3HTIPS02* 

22 1 

23D 


3,//,<15X,F6.3,F 1 3.5,Fl2.5.Fll .5 1) 

*02* 

222 

231 

37n 

NDIFF (K >=L1 

_ _ *02* 2 23 

232 

380 

IfCNSET2.EQ.O»6o To 9 5 q 

*0?* 

229 

233 


00 91,0 K = 1 , N SET 2 

f imui'il Mil 

F£Ta 

23*1 


READ(LCG1,39C)L1 ,L2 

* C2* 

226 

235 

390 

for„ati?i3) 

*02* 

227 

236 


CALL UD0303CLNCT,7«L1) 

*02* 

228 

237 


NM(K »=L1 

*02* 

.2.29 

238 


NRAD(K1-L2 

*02* 

230 

239 


READCLOgI , 400>TERao * : ,k> , COM C J, l.K) . wFR AC CJ.l .K », J=1 r Ll 1 

* r2l 

231 

290 

900 

FORMAT C FJ2.0 ,/, (2F1 2.01 1 

*02* 

232 

291 


WRITE<L0G2,4 10 »K ,L1 ,L 2, TER AD C1,K),(DMCJ,1 ,K 1 , WFRAC ( J, 1 , K ) , J=1 

,LI) *D2* 

233 

292 

910 

FORMAT ( 2X ,/ ,10X ,51HFRACTI0NAL LOSS DISTRIBUTION CURVES FOR BlADE C*o2» 

234 

293 


1LA8S»I2»I5, 16H POINTS GIVEN A T . 1 3 .1 7 H R AD I A L LOCATIONS, //,10X»52HF*02* 

235 

299 


2R ACll ON OF C0MPUTIN6 STATION LENGTH AT BLADE EXIT = ,F7 . 9 , // , 1 OX ,28 * o2* 

236 

295 


3HFR ACT 1 ON OF MERIDIONAL C HCR n *4 X , 26 hLOS S/ 1 OSS AT T RAILING EDGE,//, *02* 

237 

296 


9|15X,F11«4,2GX,F11.4 J) 

*02* 

238 

297 


IF<l2«E0. 1 )G0 TO 490 

* D2* 

23.9 

298 


DO 820 L=2,L2 

*02* 

290 

299 


CALL U003Q3ILNCT , 5«L 1 > 

*02* 

24 1 

250 


REA0CL0G1 .900ITERA0 (L,K J t (DM C J,L ,K| ,WFR ACCJ.L, Kl, J=1,L1 ) 

*02* 

242 

251 

,20 

WRITE(L0G2»43ClT£FADCL,KI,lDMtJ,L,K) ,uFRACCJ,L,K).Jrl.LI) 

*02* 


252 

930 

FORMAT (2X,/,10X, 52HFRACTI0N OF COMPUTING STATION LENGTH At BLADE E*02* 

244 

253 


1XIT r,F7.4,//,10X,28HFRACTI0N OF MERIDIONAL CHORD, 9 X,26HL0SS/L0SS *02* 

245 

259 


2AT TRAILING EDGE, //,C15x,F11.9 f 20X,Fll. 4)» 

*02* 

246 

255 

9«,0 

CONTINUE 

*02* 

247 

256 

950 

IF(NSPLIT.E0.0.AND # NREA0.E0.0»G0 TO 570 

*02* 

248 

257 


READ CL0G1 ,46 0 CDELF CJ) ,J=1,NS TRMS ) 

*02* 

249 

258 

960 

FORMAT C6F12 ,0 > 

*02* 

250 

259 


Ll = 5 

*02* 

251 

260 


IFCNSTRMS.GE.16 |L1=8 

*02* 

252 

2 61 


CALL U003o3(LNCt,L1I 

*02* 

253 

2&2 


URIT£CL0G2,47C» 

*02* 

? ' * 

263 


L 1-NSTRMS 

1 r " % 

i. • ■ 

269 


IFCNSTRMS.6T.I5 ILlrlS 

■ 02* 

2 

265 


WR1TECL0G2, 480)1 J,J=I ,L1 1 

S 02* 

257 

266 

9 80 

FORMA T C 2X, /, 10X .IOhSTREA ml INE , I 5, 1917) 

*02* 

258 

267 

9 70 

FORMAT <2X,/,10X»78HPR0P0R TIO.N S OF TOTAL FLOW BETWEEN HUB AND 

EACH *02* 

259 

268 


1STREAMLINE ARE TO BE AS FOLLOUS) 

*02* 

260 

2 69 


WRITE(L0G2,99G) (DELF (J),J=1,L1> 

*02* 

261 

270 

990 

FORMAT 1 10X,4HFL0W,7X , 1 5F7 ,9) 

*02* 

262 

271 


IFCNSTRMS.LE.15JG0 TO SDO 

*02* 

263 

Tt2 


L1=L1«I 

*02* 

264 

273 


WRITE CL0G2,460»{J,J-L1,NSTR4S» 

*02* 

265 

2 79 


WRlTECL062,9 9ol t OELF C J) , J-Ll , NS TRMS > 

*02* 

266 

275 

5 00 

IFCNREAD.EQ.C160 TO S7D 

*02* 

267 

276 


READjLOGl ,51C>CCRCJ,I),XC J,I » ,XLCJ,I ) ,1 1 1 J, I ) , J J C J . I > , J=1 , NS TR MS > ,S (j2 * 

268 

277 


11=1 ,NSTNS» 

*02* 

269 

278 

510 

F0RMATC3F12. 0,2131 

*02* 

270 

279 


CALL U00303CLNCT,S«NSTRMS» 

*02* 

271 

280 


WRITE CLOG 2, 520) 

*02* 

272 

281 

520 

F0RHAT(2x»/,l0X,32HESTIMATE0 STREAMLINE COORDINATES) 

*02* 

273 

282 


00 530 I=1,NSTNS 

SO? * 

274 

283 


IFCI.GT.DCALL U0030 3 ( LNCT, 3 ♦ NS TRMS ) 

S 02 S 

275 

289 

530 

WRITECL0G2.S90) C I , J , R CJ, I ) ,X ( J, I ) ,XL C J, I ) ,1 1 1 J ,1 ) , J JC J, I ) , J = 1 

,NSTR*02S 

276 


25 


285 


IMS) 

F OR HATt?X»/tlnX«79H STATION 


_5J 


1 l-coordin a te 

2 .15 H 


checks- I 


550 


GO TO 570 
WRIT£<L0G2,560» 


560 FORMAT (THltlGX, 33 HJ OB STOPPED 

STOP 

570 RETURN 

END 
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Source listing of the 1964 - 65 Version 5 CALCOMP subroutine, UD0312, 
appears below. 


SUBROUTI NE U00312 


* 12 * 


8 

9 

10 

11 


REAL LOSS .LAMI , LAMIPj , LAMIM 1 *12* 3 

COMMON NSTNS . NS TRMS .NMAX .NFORCE ,N B L , NCA SE .U S PL IT .NR EAO . NP UNrH . NP A p* 1 2* A 
1E.NSET1,NSET2,ISTAG,ICASE,IFAIL0,IPASS, I, IVFAXL ,IFFAIL,NMIX,NTRANS* 12 * 5 
2 , NPLOT ,ILOSS,LNCT,ITUr,IHIq,IfAtL,iT£R.| nG 1 .L 0 G 2 . L 0 G 3 .LOG A .L065 .LO* 1 2* 6 

3G6,IPRINT,NMANY,NSTPLT,NEQN * 12 * 7 

COM MON N SP E C I301 , Ny ORK 13 0 ) ,NLOS S j 30 L .N0A TA(30 1 . NTE R P 1 3n ) . NM ACH j MO > * 12* 8 

1.NL1 (3C 1,NL2(30 I .NDiMENC 30), IS! ip » , X s2 ( 3 0 I , I S 3 ( 30 1 ,NEV AL ( 30 1 ,ND 1 F* 12 * 9 

2F CM ) , NOEL 1301, N LlTE Pt 1 Q),HMI 2 1, Nft AD ( 2 1 , NC UR Vr ( 3 C 1 , N Wm CH ( 3 0 1 , NOUT 1* 12* 1 CL 

3 t 3a 1 ,N0tlT2 ( 30 1 ,N OUT 3 I 30 | » NBL ADC f 30) *12* 11 

COMMGN pHUl ,S,2 1 TRAC 1 1 1 ,5 , 2) *12c 12 


12 

13 


COMMON R«2i,30),XL(21,30l t XI21, 301, H 121 ,301, S (2 1,30 1,VN 121,301 1*12* 


13 

1* 


14 

15 

2, TeRAD(5,2) 

COMMON OAT AC 1 10 0 l.OATAH 1001 . nATA2« 1 001 .DATA 3 11 001 . 0AtA4( 1001 

*12* 
.DATA* 12* 

15 

16 

16 

17 

IS (ICC 1, OAT A6 ( 1QQ),D A TA 7(1 001 .DATA 8C 1C01 ,DATA 9 C 1 DO 1 , FL0U 1 1 0 1 .SPEED «* 12* 
2301 .SPDFAC ( 1C1,R BLOCK ( 30 1 »BD 1ST ( 301 , u BLOCK (30 1 .WUBL (301 ,XSTN( 15 a) ,* 12* 

17 

1* 

18 

19 

jRSTN ( 150 ) ,OELF ( 3 0),DELc( 1001 .DELTA ( 1 DO) , TITLE ( 1 8 1 ,D R0M2 ( 30 1 ,R IM1 1 3* 1 2* 
40 ),XIM1 ( 3 C ) . u OR K (21),L0SS(2i 1 ,T ANEPS (21 1.XK21 ),WV( 21),0£LU(21),i A*12* 

19 

2Q 

20 

21 

5mI( 211.LAMIM1(21),LAMIP1(21).PHT(21),CR(211,GAMA(211,SPPG(21) 
6(21 ),HKFEP(211,SKrfrP(21) ,Vk<EEP (21) ,DELH( 30).DELT(3C1 

,CPPG* 12* 
* 12* 

21 

22 

22 

23 

COMMON VlSK .SHAPE .SCLFAC , E J , G .TOLNCE , XS C* LE ,PS C AL E . PLOU ,R L 0U , 
1. RCONST .FM2.HM1N ,C1 ,PT ,CONTR .CONMX 

XhmAXS 12* 
*12* 

23 

29 

24 

DIMENSION PSTAT ( 321 ,XX(32 1 

*12* 

25 

25 

X M A X -X ( 1 .NSTNS) 

*12* 

26 

26 

27 

XMI N-X ( 1,1) 

□0 100 J-2, nSTRmS 

*12* 
* 12* 

27 

2.8 

28 

IF ( X (J, 1).LT.XHIN)XMIN=X ( J,1 » 

*12* 

29 

29 

IF ( x (J .NSTNS l.GT.XMAX )XMAX=X ( J, NSTNS 1 

*12* 

30 

30 
3 1 

100 CONTINUE 

XMIN=FLOaT(IFIX(XMxN)) 

*12* 

MOD. 

31 

-JGU 

32 

IF( (XMIN.GE.O. 0) .ANO. (XMIN.LE . 5 . 0 1 1 X M IN =0 ,Q 

MOD. 

-JGu 

33 

ALEN=( xMAx-XMIN 1 /XSCALE 

MOD. 

-JGU 

34 

XLEN-2C.0 

MOD. 

-JGU 

35 

IF(ALEN.LE.10.aiXLEN=10.0 

MOO . 

- J V .. 

36 

IF(NPLOT.E0.2)GO TO 134 

* 125 

92 

37 

CALL PLOTID 

MOD. 

-JGU 

38 

call plotig.o.-12.o,-3) 

MOD. 

-JGU 

39 

CALL PL0T(Q. O.o. 5,-31 

MOD. 

-JGU 


MO 

41 


CALL AXIS(i3.a,c.0,l6HAXlAL C 0 OR 0 I NA T E ,- 1 6 , XLeN . 0 • 0 , XM IN , X S CALE , 1 0 •« OD . 
10 I MOD . 


-JGU 

-JGU 


42 

43 

44 


call axisic.o.o.o.ishstatic pre ssure ,15 .io.o.ro.o.plou ,pscale ,io.omod. 

1 1 M OD . 


J = 1 


* 125 


-JGU 

-JGU 

47 


45 


K -1 

*12* 

48 

46 


7PEN-3.4 

MOD . 

-JGU 

47 


X PE N = (X MAX/ XSCALE 1-3.0 

MOD. 

-JGU 

48 


S I2E=0. 14 

MOD. 

-JGU 

49 

110 

00 120 1=1, NSTNS 

* 12* 

49 

50 


HS=h«J.I)-(VU(J, I)**2*VM( J.r 1**2 1/(2 ,0 *G*EJj 

*12* 

50 

51 


IFIHS ,LT.HMIN)HS=HMIN 

*12* 

51 

52 


PSTAT (I )=UDG4(HS ,S l J , 1 1 1 /SCL F AC **2 

* 12* 

52 

53 

120 

XX(I)=X(J,I) 

* 125 

53 

54 


CALL LINE (XX ,PST AT, NSTNS.1.1 , A, XMIN, XSC ALE.PLOU.PSCALE 1 

MOD. 

-JGU 

55 


CALL SYMBOL (XPEN , YPEN*C. 0 7 , S I Zf ,K ,0 . 0 ,- 1 1 

MOD. 

-JGU 

56 


IF( J.E0.NSTRMS1CALL S YM B 0 L ( X P EN *0 .1 4 , TP EN ,S IZE , 1 7H= > TIP 

STREAMLINMOD. 

-JGU 
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hod.-jgu 


1E,0.Q,17j 

IF( J.EQ .NSTRHS160 TO 133 

K -K ♦ 1 

IF( J.EQ.IHIDIJrNSTRHS 

IF ( J.EO .NSTRHS 1C ALL S YHB OH X PEN -*0.111 »YP EN ,SIZe ,17h= 
IE ,0.0,171 


IF* J.fO.l 

IFtJ.EO.lHIO ICALL S VHBOL < XPEN +0 . t 4 £ V P , ,N 

10. □ ,17 ) 

YPeN—T PEN-0. 25 

GO TO HO 

130 CALL PLOT(25.n,-17.Q,- 
IFtNPLOT.EQ.llGO TO 18C 

1 3 *« CALL PLOT fp 

CALL PLOTtO. 0,-12. 0,-31 

CALL PLOTtO. G,Q. 5,-31 ■ 

CALL AXlS(0.0,a.o,i6HAXIAL C 0 CR D I NA TE ,- 16 ,XLeN , 0 .0 , 
ID » 


C ALL AXIS ( 0.0, 0.0,6 HR AOI US, 6 , 10 »D»90 .0, RLOw.XSCALE, 

DO 150 J=i, NST RHS 

DO 140 Izl.NSTNS 
X X ( T 1-X (J,I1 


14o PSTATtI)-R(J,Il 

150 CALL LINE (XX .PSTAT.NSTNS, t , n # 11 , X MI N , XSCALE.RLOU.XS 


DO 170 I=1,NSTNS 

p O 160 .|=1 , NSTRhS 

PSTATtJJ -R ( J , 1 1 

160 XXtJl=X(J,Il 

170 CALL LInE(XX,p st At,NSTRHS,1,0,11,XHIN,XSCALE,RL0U,X 
CALL PLOT (25.0,-12.0,-31 


180 CALL ENDCC 

RETURN 

END 


tl2t 

*12* 

> MiD STREAHLlNHOD.- 
hod.- 


* 12 * 


*125 62 

HO O.-J SU 

HOD.-JGU 

HOD.-JGU 

XHiN.XSCAlE.Iq. HOD.-JG u 


10.0 | 


SCALE1 


H OD . -J 

GU 

MOD.-J 

*12* 

GU 

67 

512S 

mt 

68 

6SL 

*12* 
MflO .“J 

70 

GU 

*12* 
5 12* 

76 

77 

*12* 

*12* 

70 

79 

MOD.-J 

HOO.-J 

GU 

GU 

MOD.-J 

GU 
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Axial coordinate 

Figure 2. - Definition of meridional direction. 



Figure 4. - Fan module axial station locations. 







O Tgfl'i 



Percent span from hub 


Figure 9. - Temperature profiles at axial stations. 
Ti ■ 442. 82K (797. 08R). 


















NORMALIZATION CONSTANT -133.47 M/SEC 



(b) Station a 
Figure 11 - Concluded. 






NCRMGUlZfiTION CONSTANT -33. j; 


C -> INOUSTFIT PRCCRfiK OUTPUT 
& -* MODIFIED PROGRAM OUTPUT 



( } 1 1 1 1 - 

3 40.00 50.00 60.00 70.00 90.00 90 . 

PERCENT SPAN PROM HUB 

(b) Station 8. 

Figure 14 - Concluded. 








Figure 16. - Streamline contraction through fan module (program generated CALCOMP plot). One unit of length 
on the axes equals 2.54 cm. (LO inch). 



RXIflL COORDINATE 


Figure 17. - Static pressure distribution in fan module (program generated CALCOMP plot). One unit of pressure 
equals 6. 895 kPa (LO psia). One unit of length equals 2 54 cm. (LO inch). 
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